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<^ Abstract 

§^ First observations of the 6° ^ iK2S)ti, 6° ^ •ilj(2S)7r+7t- and 6° ^ •iI;(2S)7T+7r- de- 

^ cays are made using a dataset corresponding to an integrated luminosity of 1.0 fb~^ 

coUected by the LHCb experiment in proton-proton coUisions at a centre-of-mass 
energy of ^/s = 7TeV. The ratios of the branching fractions of each of the ij^l^S) 
modes with respect to the corresponding J/t|) decays are 



X 



0.83 ± 0.14 (stat) ± 0.12 (syst) ± 0.02 (B), 



> 

J^i g(B!^^r|;(2S)Ti) 

^ S(BO ^ J/rl;r|) 

^ B(B^ — )• tbf2S')7t+7T~) 

^TF^n — _^ = 0.56 ± 0.07 (stat) ± 0.05 (syst) ± 0.01 (S), 
0(0*^ — J/i|;7r+7T ) 

^ = 0-34 ± 0.04 (stat) ± 0.03 (syst) ± 0.01 (S), 



i3(B0 ^ J/i|j7r+7T-) 

where the third uncertainty corresponds to the uncertainties of the dilepton branch- 
ing fractions of the J/i|j and ^i'(2S) meson decays. 
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1 Introduction 



Decays of B mesons containing a charmonium resonance, J/i|) or i|;(2S), in the final state 
play a crucial role in the study of CP violation and in the precise measurement of neutral 
B meson mixing parameters. 

The Bg — )■ J/il^ri decay was observed by the Belle collaboration and the branching 
fraction was measured to be -B(B° J/il^rj) = (5.10±0.50±0.25 +^;^^) x 10"^ [l], where the 
first uncertainty is statistical, the second systematic and the third due to the uncertainty 
in the number of produced B^B^ pairs. This decay has also recently been reported by 
LHCb, including the decay B° ^ J/ipri' |2j. 

The B°g^ — )■ J/i|; Tt+Tt" decays, where B°g^ denotes a B° or B^ meson, have been studied 
previously and the W^tC final states are found to comprise the decay products of the 
p°(770) and f2(1270) mesons in case of B° decays and of fo(980) and fo(1370) mesons in 
case of Bg decays [3]-[5]. The Bg modes have been used to measure mixing-induced CP 
violation j6,7j. The decays B^ — > \l;(2S)ri and B^^^ — > \l;(2S)7r"'"7t~ have not previously 
been studiecL 

The relative branching fractions of B'^ and B^ mesons into final states containing J/ij) 
and i|'(2S) mesons have been studied by several experiments (CDF |8|[9], DO |lO] and 
LHCb |ll]). In this paper, measurements of the branching fraction ratios of B9 mesons 
decaying to ij»(2S)X° and J/ij^X" are reported, where X° denotes either an r| meson or a 
Tt+Tt" system. Charge conjugate decays are implicitly included. The analysis presented 
here is based on a data sample corresponding to an integrated luminosity of 1.0 fb~^ 
collected with the LHCb detector during 2011 in pp collisions at a centre-of-mass energy 
of v^= 7TeV. 



2 LHCb detector 



The LHCb detector |12] is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < 77 < 5, designed for the study of particles containing b or c quarks. The 
detector includes a high precision tracking system consisting of a silicon-strip vertex de- 
tector surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4Tm, and three stations 
of silicon-strip detectors and straw drift tubes placed downstream. The combined track- 
ing system has momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% 
at 100 GeV/ c, and impact parameter resolution of 20 (xm for tracks with high transverse 
momentum (px)- Charged hadrons are identified using two ring-imaging Cherenkov de- 
tectors. Photon, electron and hadron candidates are identified by a calorimeter system 
consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter 
and a hadronic calorimeter. Muons are identified by a system composed of alternating 
layers of iron and multiwire proportional chambers. 



The trigger 13 consists of a hardware stage, based on information from the calorimeter 



and muon systems, followed by a software stage where a full event reconstruction is 



1 



applied. Candidate events are first required to pass a hardware trigger which selects 
muons with a transverse momentum, px > 1-48 GeV/c. In the subsequent software trigger, 
at least one of the final state particles is required to have both px > 0.8 GeV/c and impact 
parameter > 100 |J.m with respect to all of the primary pp interaction vertices (PVs) in 
the event. Finally, two or more of the final state particles are required to form a vertex 
which is significantly displaced from the PVs. 

For the simulation, pp collisions are generated using Pythia 6.4 ^lAj with a specific 
LHCb configuration [15j. Decays of hadronic particles are described by EvtGen |16| 



in which final state radiation is generated using Photos 17 . The interaction of the 



generated particles with the detector and its response are implemented using the Geant4 



toolkit 18,19 as described in Ref. 120 



3 Event selection 

The decays B^^^ — > iprj and B^^,^ — t- i|>7r+7r^, where ij) denotes J/ip or 11^(28), are recon- 
structed using i|) — 7- and r[ — )• yy decay modes. Pairs of oppositely-charged tracks 
identified as muons, each having px > 0.55 GeV/c and originating from a common vertex, 
are combined to form \|; — )■ jj.+ ix" candidates. Track quality is ensured by requiring the 

per number of degrees of freedom (x^/ndf) provided by the track fit to be less than 5. 
Well identified muons are selected by requiring that the difference in logarithms of the 
global likelihood of the muon hypothesis, Alog£^h (2l], provided by the particle identifi- 
cation detectors, with respect to the hadron hypothesis is larger than zero. The fit of the 
common two-prong vertex is required to satisfy x^/^df < 20. The vertex is deemed to 
be well separated from the reconstructed primary vertex of the proton-proton interaction 
by requiring the decay length significance to be larger than three. Finally, the invariant 
mass of the dimuon combination is required to be between 3.020 and 3.135 GeV/c^ for J/ij) 
candidates and between 3.597 and 3.730 GeV/c^ for ■i['(2S) candidates. These correspond 
to [-5cr; 3cr] intervals around the nominal masses to accomodate QED radiation. 

The pions are required to have px > 0.25 GeV/c and an impact parameter x^, defined 
as the difference between the of the PV formed with and without the considered track, 
larger than 9. When more that one PV is reconstructed, the smallest value of impact 
parameter is chosen. In addition, to suppress contamination from kaons, the difference 



between the logarithms of likelihoods of the pion and kaon hypotheses, A log C^k 22 
provided by the RICH detectors, has to be larger than zero. 

Photons are selected from neutral clusters in the electromagnetic calorimeter with 
transverse energy in excess of 0.4 GeV. The r\ yy candidates are reconstructed as dipho- 



ton combinations with an invariant mass within ±70MeV/c of the r| mass 23 . To sup- 
press the large combinatorial background from the decays of neutral pions, photons that 
form a 71*^ — yy candidate with invariant mass within ±25 MeV/ of the tjP mass are not 
used to reconstruct r[ — )■ yy candidates. 

The B^g^ candidates are formed from i|)X° combinations. In the iprj case an addi- 
tional requirement Pt{t]) > 2.5GeV/c is applied to reduce combinatorial background. To 
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improve the invariant mass resolution a kinematic fit [24j is performed. In tfiis fit, con- 



straints are applied on the known masses 23 of intermediate resonances, and it is also 



required that the candidate's momentum vector points to the associated primary vertex. 
The x^/iidf for this fit is required to be less than 5. Finally, the decay time, ct, of the 
B°g^ candidate, calculated with respect to the primary vertex, is required to be in excess 
of 150 (xm. 



4 Observation of the i|)(2S)r| decay 

The invariant mass distributions of the selected ^\)r\ candidates are shown in Fig. [ij The 
Bg — i- '\\)r\ signal yields are estimated by performing unbinned extended maximum like- 
lihood fits. The Eg signal is modelled by a Gaussian distribution and the background 
by an exponential function. In the J/i|>ri case a possible contribution from the corre- 
sponding B° decays is included in the fit model as an additional Gaussian component. 
The resolutions of the two Gaussian functions are set to be the same and the differ- 
ence of their central values is fixed to the known difference between the and the B° 
The contribution from the decay B'' — )■ \[)(2S)ri is not considered in the base- 



masses 



23 



line fit model. The mass resolution of the B^ — )■ \[)(2S)ri decay mode is fixed to the value 



^|.(2S)rl 
DATA 



^DATA ^ '^MC 



I ^uc ' where o"data and umc are the widths of the corresponding 



channel in data and simulation, respectively. 
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Figure 1 
function 
show the 



: Mass distributions of (a) B'^ 



J/\1;ti and (b) B' 



(s) 



\[>(2S)ri candidates. The total fit 
(solid black) and the combinatorial background (dashed) are shown. The solid red lines 
signal Bg contribution and the red dot dashed line corresponds to the B'^ contribution. 



The fit results are summarised in Table [T) In all cases the positions of the signal peaks 
are consistent with the nominal B^ mass 23 and the resolutions are in agreement with 



the expectations from simulation. The measured yield of B° — )■ J/v|jri is 144 ± 41 events 
(uncertainty is statistical only), which is consistent with the expected value based on the 
measured branching fraction of this decay i25j . The statistical significance in each fit is 
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Table 1: Fitted values of signal events (A^b); signal peak position (Mb) and resolution (cr^). 
The quoted uncertainties are statistical only. 



Mode A^B 

B^^TlAini 863 ± 52 
1^(28)11 76 ± 12 



Mb ctb 

[MeV/c^] [MeV/c^] 

5370.9 ±2.3 33.7 ±2.3 
5373.4 ± 5.0 26.6 fixed 



determined as iS = y — 21n -^^^ , where £s+b and Cb denote the likelihood of the signal 

plus background hypothesis and the background only hypothesis, respectively. Taking 
into account the systematic uncertainty related to the fit function, which is discussed in 
detail in Sect. |6| the significance of the — )■ i|;(2S)ri signal is 6.2cr. 

To demonstrate that the signal originates from — > \|j(2S)ri decays the sPlot tech- 
nique [26j has been used to separate the signal and the background. Using the |j."'"|j.~yy 
invariant mass distribution as the discriminating variable, the distributions for the in- 
variant masses of the intermediate resonances r\ — t- yy and 'i|'(2S) — t- M-^M-^ have been 
obtained. In this procedure, the invariant mass window for each corresponding resonance 
is released and the mass constraint is removed. The resulting invariant mass distributions 
for yy and \Ji^[i~ from B^ — > \l;(2S)ri candidates are shown in Fig. |2} Clear signals are 
seen in both r\ — )• yy and i|'(2S) — > \L^\^.~ decays. The distributions are described by the 
sum of a Gaussian function and a constant. The fit shows that the constant is consistent 
with zero, as expected. 



LHCb 




0.46 0.48 0.5 0.52 0.54 0.56 0.58 0.6 0.62 0.64 

M(yy) [ GeV/c2] 




3.62 3.64 3.66 3.68 3.7 3.72 3.74 



M(^+^l-) [GeV/c2 



Figure 2: Background subtracted (a) yy and (b) ]x^\i mass distributions in 6° — ^ T[)(2S)ri 
decays. In both cases the blue line is the result of the fit described in the text. 
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Table 2: Fitted values of signal events (A^b); signal peak position (Mb) and resolution (cr^). 
The quoted uncertainties are statistical only. 

l^/r A AT CTb 

^"^^ [MeV/c^] [MeV/c^] 

B° ^ J/il)7t+7t- 2801 ±85 5281.1 ±0.3 8.2 ± 0.3 

B° J/il)7t+7t- 4096 ±86 5368.4 ±0.2 8.7 ±0.2 

B° ^ •iK2S)7r+7t- 202 ±23 5280.3 ± 1.0 8.4 ±1.1 

B° ^ •iIj(2S)7r+7t- 178 ±22 5366.3 ±1.2 9.1 ± 1.4 

5 Observation of the B^^^ i()(2S)7r+7r~ decays 

The invariant mass distributions for the B^^^ — )■ T|)7r+7t^ candidates are shown in Fig. 3j 
The narrow signals correspond to the B'' — )■ i|)7r"'"7t~ and B^ — )■ i|)7r+7r~ decays. The peak 
at lower mass corresponds to a reflection from B° — )■ \|jK*°(— )■ K"^7r~) decays where the 
kaon is misidentifled as a pion. The contribution from B^ — )■ i|)K*° decays [27| is negligible. 

The invariant mass distributions are fltted with two Gaussian functions to describe 
the two signals, an asymmetric Gaussian function with different width for the two sides 
to represent the reflection from B" — )• \[)K*° decays and an exponential function for the 
background. The flt results are summarised in Table [2j The statistical significances of 
the signals are found to be larger than 9 standard deviations. 

For the B^^^ — )■ J/i|)7r+7t^ decays, the tt+tt^ mass shapes have been studied in detail us- 
ing a partial wave analysis in Refs. |4|5]. The main contributions are B° — )■ J/i|; p°(770) and 




M(J/i|;7r+7t-) [GeV/c2] M(tK2S)7T+7T-) [GeV/c^] 



Figure 3: Mass distributions of (a) B^^^ — ;> J/\1j7T+7T and (b) 6°^^ — )- ^\){2S)n~^7X candidates. 
The total fit function (solid black) and the combinatorial background (dashed) are shown. The 
solid red lines show the signal contribution and the red dot dashed lines correspond to the B*^ 
contributions. The reflections from misidentified B*^ — >■ tJjK*'^, K**^ — )• K^7r~ decays are shown 
with dotted blue lines. 
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M(7t+7T-) [GeV/c2] M(7t+7T-) [GeV/c^] 



Figure 4: Background subtracted tt+tt" mass distribution in (a) B° — )• •ij;(2S)7T^7r~ and 
(b) Eg — )• r|j(2S)7t+7t^ (black points). The red filled area shows the expected signal spectrum 
for the iK2S) channel derived from the measured spectrum of the J/4) channel (the fit has one 
parameter — the normalisation). The width of the band corresponds to the uncertainties of the 
distribution from the J/4» channel. In case of B° •il;(2S)7T+7r-, the blue vertical filled area 
shows the Kg region that is excluded from the fit. 



Bg — )■ J/il;fo(980). However, due to the limited number of signal events, the same method 
cannot be used for the B^^,^ — )■ i|)(2S)7t^7r^ decays. The sPlot technique is used in order to 
study the dipion mass distribution in those decays. With the \[>(2S)7r"'"7r~ invariant mass 
as the discriminating variable, the tt+tt" invariant mass spectra from B°j,j — > \[)(2S)7r+7t~ 
decays are obtained (see Fig. 

To check that the background subtracted Tt^n^ distributions have similar shapes in 
both channels, the distribution obtained from the i|»(2S)7r+7r^ decay is fitted with the 
distribution obtained from the J/i[)7r"^7r~ channel, corrected by the ratio of phase-space 
factors and by the ratio of the efficiencies which depends on the dipion invariant mass. 
The p- value for the fit is 30% for B° — )■ ij;7t"''7r~ and 7% for Bg — )■ '\\)n'^7i~ , respectively. 
As seen in Fig. g B° -> •iIj(2S)p°(770) and B° ^ iIj(2S)fo(980) decays are the main 
contributions to B^^^ \[)(2S)7t+7t^ decays. Detailed amplitude analyses of the resonance 
structures in B^^^ — )■ \[)(2S)7r+7t^ decays, similar to Refs. j4 5j, will be possible with a larger 
dataset. This will allow the possible excess of events in the region M(7t^7r^) > 1.4GeV/c^ 
to be investigated. 

The narrow peak around 0.5GeV/c^ in Fig. |4]^a) is dominated by Kg — )■ W^n' from 
B*^ — )• J/4) Kg decays. The contributions from Kg decays are taken into account by the fit 
function described in Ref. [2]. The resulting yields are 129 ± 26 in the J/4) channel and 
11 ± 6 in the ■4'(2S) channel. In the calculation of the final ratio of branching fractions, 
the number of Kg events is subtracted from the corresponding B° — )■ xI^tt+tt" yields. The 



yield from B^ — )■ ij^Kg decays is negligible 28 
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6 Efficiencies and systematic uncertainties 



The ratios of branching fractions are calculated using the formula 
g(B ^ -i|;(2S)X°) ^ iV^(2S)xo ^ 6j/^xo g(J/i|; ^ 

where is the number of signal events, and e is the product of the geometrical acceptance, 
the detection, reconstruction, selection and trigger efficiencies. The efficiency ratios are 
estimated using simulation for all six decay modes. 

The efficiency ratios are 1.22 ± 0.01, 1.03 ± 0.01 and 1.02 ± 0.01 for the B° 
B*^ — ^^T^^^l~ and B^ — > '\]^n^TC channels, respectively (uncertainties are statistical only). 
Since the selection criteria for the decays with J/i[) and iK2S) are identical, the ratio of 
efficiencies is expected to be close to unity. The deviation of the overall efficiency ratio 
from unity in case of B^ — )■ ijiri is due to the difference between the spectra of the 
selected J/ij) and 'i|)(2S) mesons, when the PtIti) > 2.5GeV/c requirement is applied. For 
the B^g^ — )■ i|)7T"'"7T~ channels this effect is small since no explicit requirement is applied 
on the dipion system. 

Most systematic uncertainties cancel in the ratio of branching fractions, in partic- 
ular, those related to the muon and reconstruction and identification. Systematic 
uncertainties related to the fit model are estimated using a number of alternative mod- 
els for the description of the invariant mass distributions. For the B^ -> ij;ri decays 
the tested alternatives are a fit model including a B*^ signal component (with the ratio 
iV(B° — )■ i|)r|)/iV(Bg — )■ i|;ri) fixed from the J/ij) channel), a fit model with a linear function 
for the background description, fits with signal widths fixed or not fixed to those obtained 
in simulation, a fit with the difference between the fitted B° and B? masses allowed to 



vary within a ila interval around the nominal value 23 , and a fit model with Student's 
t-distributions for the signals. For each alternative fit model the ratio of event yields is 
calculated and the systematic uncertainty is then determined as the maximum deviation 
of this ratio from the ratio obtained with the baseline model. For B^^^ — ipTt+Tt" decays 
the tested alternatives include a fit with a first or second order polynomial for the back- 
ground description, a model with a symmetric Gaussian distribution for the reflection and 
a model with the difference of the mean values of the two Gaussian functions flxed to the 
known mass difference between the B^ and the B° mesons [23j. The maximum deviation 
observed in the ratio of yields in the ■i['(2S) and J/i|» modes is taken as the systematic 
uncertainty. The obtained uncertainties are 8.0% for the B^ — )■ ipr] channel, 1.0% for the 
B° — 7- ^\>n'^n~ channel and 1.6% for the B^ — )■ ipTr+Tr" channel. 

The selection efficiency for the dipion system has a dependence on the dipion invariant 
mass. The ratios of efficiencies vary over the entire ti^tC mass range by approximately 
40% and 24% for B° — )■ ipTt+Tt^ and B^ — )■ i|;7r+7r~ channels, respectively. The systematic 
uncertainties related to the different dependence of the efficiency as a function of the 
dipion invariant mass for J/ij; and 1K2S) channels are evaluated using the decay models 
from Ref. js] for Bg and Refs. [2|[4] for B*^ decays. The systematic uncertainties on the 
branching fraction ratios are 2% for both channels. 
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The most important source of uncertainty arises from potential disagreement between 
data and simulation in the estimation of efficiencies. This source of uncertainty is studied 
by varying the selection criteria in ranges corresponding to approximately 15% change in 
the signal yields. The agreement is estimated by comparing the efficiency corrected ratio 
of yields with these variations. The resulting uncertainties are found to be 11.5% in the 
Bg — )■ \|;ri channel and 8% in the B^^^ — )■ '\\)n~^n~ channel. 

The geometrical acceptance is calculated separately for different magnet polarities. 
The observed difference in the efficiency ratios is taken as an estimate of the systematic 
uncertainty and is 1.1% for the B° — t- \|)7t+7t~ channel and negligible for the other channels. 

The trigger is highly efficient in selecting B meson decays with two muons in the final 
state. For this analysis the dimuon pair is required to trigger the event. Differences in 
the trigger efficiency between data and simulation are studied in the data using events 



that were triggered independently of the dimuon pair 11 . Based on these studies, an 
uncertainty of 1.1% is assigned. A summary of all systematic uncertainties is presented 
in Table [3l 



Table 3: Relative systematic uncertainties (in %) of the relative branching fractions. 



Source 




B° ipTt+Tr 


Bg — > i|;7r+7r 


Fit model 


8.0 


1.0 


1.6 


Mass dependence of efficiencies 




2.0 


2.0 


Efficiencies from simulation 


11.5 


8.0 


8.0 


Acceptance 


< 0.5 


1.1 


< 0.5 


Trigger 


1.1 


1.1 


1.1 



Sum in quadrature 14.1 8.5 8.5 



7 Results 

With data corresponding to an integrated luminosity of 1.0 fb~^, collected in 2011 with the 
LHCb detector, the first observations of the B^ — )■ ij)(2S)ri and B^^^ — )■ i|)(2S)7t"'"7t~ decays 
have been made. The relative rates of B^^^ meson decays into final states containing 
J/4> and ij»(2S) mesons are measured for those decay modes. Since the dielectron branching 
fractions of ij) mesons are measured more precisely than those of the dimuon decay modes, 
invoking lepton universality, the ratio ^^^[23^^^,+^,-) = g^|f?sHe+e-) = 7.69 ± 0.19 



23 



IS 
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used. The results are combined using Eq. ([T]), to give 




where the first uncertainty is statistical, the second systematic and the third from the 



states. The branching fraction ratios measured here correspond to the time integrated 
quantities. For the B° 3/\\){\\){2S))n+n- channel the measured ratio excludes the 
Kg — > 71^ 7t^ contibution. The dominant contributions to the B^^,^ — )• i|;(2S)7t"'"7t~ decays 
are found to be from B° ^ iIj(2S)p°(770) and B° ^ •iIj(2S)fo(980) decays. 

These results are compatible with the measured range of relative branching fractions 
of B decays to -^{28) and JAp mesons. The B° Ti'(2S)ri and B° \|;(2S)7r+7r" decays 
are particularly interesting since, with more data becoming available, they can be used 
to measure CP violation in Bg mixing. 
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